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1. INTRODUCTION 

A crude estimate of the gravitational luminosity of an object of mass M, mean 
radius R and internal velocities of order V can be derived from the quadrupole 
formula 

where Rs := 2GM /(? is the Schwarzschild radius associated with the mass M 
and s is some asymmetry factor: s = for a spherically symmetric object and 
s ^ 1 for an object whose shape is far from that of a sphere. According to 
formula (|l|), the astrophysical objects for which s ^ 1, R ^ Rg and V ^ c 
may radiate a fantastic power in the form of gravitational waves: L ^ /G = 
3.6 X 10^^ W, which amounts to 10^^ times the luminosity of the Sun in the 
electromagnetic domain! 

A neutron star has a radius quite close to its Schwarzschild radius: R ~ 
1.5 — 3i?s and its rotation velocity may reach V ^ c/2 at the equator, so that 
they are a priori valuable candidates for strong gravitational emission. The cru- 
cial parameter to be investigated is the asymmetry factor s. It is well known 
that a uniformly rotating body, perfectly symmetric with respect to its rotation 
axis does not emit any gravitational wave (s = 0). Thus in order to radiate 
gravitationally a neutron star must deviate from axisymmetry. P. Haensel's 
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lecture |2| investigates the deviation from axisymmetry resulting from irregu- 
larities ( "mountains" ) in the solid crust or from the neutron star precession. 
In the present lecture, we investigate two other mechanisms which generate a 
deviation from axisymmetry: (i) the spontaneous symmetry breaking result- 
ing from the development of a triaxial instability in a rapidly rotating neutron 
star (§ ^ and (ii) the distortion induced by the internal magnetic field of the 
neutron star (§ |3|). 

2. SPONTANEOUS SYMMETRY BREAKING 

A rotating neutron star can spontaneously break its axial symmetry if the ratio 
of the rotational kinetic energy T to the absolute value of the gravitational 
potential energy, \W\, exceeds some critical value. This may occur in two 
different astrophysical situations: a just born neutron star resulting from a 
supernova may accrete matter that has not been ejected by the shock wave, 
thereby increasing its kinetic energy; alternatively, an old neutron star in a 
close binary system, accreting matter steadily from its companion, may be 
spun up until the ratio T/IVFj is high enough to allow the symmetry breaking. 

When the critical threshold r/|V7| is reached, two kinds of instabilities may 
drive the star into the non-axisymmetric state: 

1. the Chandrasekhar- Friedman- Schutz instability (hereafter CFS instabil- 

lH I Hi ; lH driven by the gravitational radiation reaction. 

2. the viscosity driven instability 1^. 

The present lecture puts the accent on the viscosity driven instability. 

2.1. Review of classical results about Maclaurin/Jacobi ellipsoids 

Let us recall some classical results from the theory of rotating Newtonian 
homogeneous bodies. It is well known that a self-gravitating incompressible 
fluid rotating rigidly at some moderate velocity takes the shape of an ax- 
isymmetric ellipsoid: the so-called Maclaurin spheroid. At the critical point 
T/|T/F| = 0.1375 in the Maclaurin sequence, two families of triaxial ellipsoids 
branch off: the Jacobi ellipsoids and the Dedekind ellipsoids. The former are 
triaxial ellipsoids rotating rigidly about their smallest axis in an inertial frame, 
whereas the latter have a fixed triaxial figure in an inertial frame, with some 
internal fluid circulation at constant vorticity (see ref. ^ or Q for a review of 
these classical results). The Maclaurin spheroids are dynamically unstable for 
T/\W\ > 0.2738. Thus the Jacobi/Dedekind bifurcation point T/\W\ = 0.1375 
is dynamically stable. However, in presence of some dissipative mechanism 
such as viscosity or gravitational radiation (CFS instability) that breaks the 
circulation or angular momentum conservation, the bifurcation point becomes 
secularly unstable against the I — 2, m — 2 "bar" mode. Note also that a 
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non-dissipative mechanism such as the one due to a magnetic field with a com- 
ponent parahel to the rotation axis breaks the circulation conservation Q and 
may generate a spontaneous symmetry breaking. If one takes into account 
only the viscosity, the growth of the bar mode leads to the deformation of the 
Maclaurin spheroid along a sequence of figures close to some Riemann S el- 
lipsoids|(^ and whose final state is a Jacobi ellipsoid Q . On the opposite, if 
the gravitational radiation reaction is taken into account but not the viscosity, 
the Maclaurin spheroid evolves close to another Riemann S sequence towards 
a Dedekind ellipsoid . 

As we shall see in the next section, this symmetry breaking is a particular 
case of a more general phenomenon . 

2.2. Spontaneous breaking of symmetry: a general phenomenon 

2.2.1. A toy model 

Let us start with a simple example. Consider the toy model of Fig. ^ : two 
heavy points of mass m, rigidly connected and sliding on a circular guide of 
radius R, of negligible mass and rotating with the angular velocity fl. It is 
straightforward to write down the expression for the potential energy of this 
system and to find the equilibrium configuration in a frame rotating with the 
guide. With the notations of Fig. |l|, the potential (potential energy per unit 
mass) writes: 

U{0) = -(i?f^)2(sin2a + cos2a sin^ 6) - 2gRcos a cosy , 2) 

g being the gravity acceleration. In spite of the invariance of the potential U 
with respect the transformation 6 ^—^ —9, two equilibium configurations exist 
if the angular velocity is larger than the value Ocrit = [ g cos a/ (i? cos 2a) ]^/^ : 
9 = and cos6' = gcosa/ {RH,^ cos 2a). 

This behaviour is typical of the phenomenon of spontaneous symmetry break- 
ing. More generally, a dynamical system is said to break spontaneously its 
symmetry if there exist solutions with lower symmetries than those of the 
Lagrangian from which the equations of motion have been derived. Other 
examples of spontaneous symmetry breaking are the phase transitions (ferro- 
magnetism at the Curie point, melting of solids), the famous walking stick in 
Charlie Chapling movies (the stick is axially symmetric but when the pressure 
of the hand on its top is larger than a critical value, the stick bends), the 
buckling collapse of metallic structures and so on. 

Coming back to the toy model, it is easy to see that, when two equilibrium 
configurations exist, the one is unstable and the other one is stable. For what 

(^) The Riemann 5 family is formed by homogeneous bodies whose fluid motion can 
be decomposed into a rigid rotation about a principal axis and a uniform circulation 
whose vorticity is parallel to the rotation vector. Maclaurin, Jacobi and Dedekind 
ellipsoids are all special cases of Riemann S ellipsoids (for more details, cf. Chap. 7 
of reL §] or Sect. 5 of reL @). 
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Fig. 1. — Toy model illustrating the phenomenon of spontaneous symmetry breaking: 
the configuration depicted with dashed lines violates the symmetry with respect to 
the rotation ajcis and is a stable equilibrium configuration if O > ficrit- 



follows it is more convenient to take the angular momentum as a free parameter 
instead of the angular velocity. Let us suppose that some viscous force acts 
on the two mass points sliding on the circular guide. If at time t = 0, the 
system is in the equilibrium configuration 9 = and its angular momentum L 
is larger than the critical value Lcrit corresponding to ficrit, then the system 
will move to the second equilibrium configuration. In this example, we have 
all the main ingredients of the spontaneous symmetry breaking mechanism: a 
parameter describing the equilibrium configuration of the system, a dynamical 
quantity that is not conserved (the kinetic energy) and a total energy that must 
be at its minimum. In order to understand the existence of two equilibrium 
configurations, it is worth writing the total energy (per unit mass) with respect 
to the inertial frame in terms of the angular momentum L: E = 1/2 L^// + 
Ugrav, where I {9) is the moment of inertia of the system with respect to the 
rotation axis and f/grav is the gravitational potential: t/grav = —2gRcosacos9. 
It is clear that, at fixed L, L^// decreases when I increases, i.e. when 9 
increases. On the contrary, Ugrav is an increasing function of 9. When L is 
large enough, the first term prevails and E decreases when 9 increases. 

In order to see the link existing between the symmetry breaking mecha- 
nism of the above toy model and the second order phase transition theory, 
let us briefly recall the Landau theory of these phenomena (see Chap. 14 of 
13 ). In a second order phase transition, a thermodynamical system makes 
a continuous transition from a state with a given symmetry to a state with a 
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lower symmetry. The leading idea in the Landau theory is the expansion of 
the thermodynamical potential <i>(P, T) in a series of a parameter rj (the order 
parameter) which measures how far the thermodynamical system is from the 
state of higher symmetry: 

$(P, T, Tj) = <i>o{P, T) + aTj + Ari^ + Bri^ + Crf + ■ ■ ■ (3) 

From invariance considerations, a must be identically equal to zero, B must 
also be equal to zero in order to obtain a stable system. On the contrary, C 
must be positive. A must be larger than zero in the symmetric phase because 
the system is stable. On the contrary, in the less symmetric phase, values of 
77 different from zero correspond to a stable state. That is possible only if 
A < Q, therefore A must vanish at the critical temperature: A = a{T — Tc). In 
conclusion we can write: 

$(P, T, 77) = T) + a{T ~ T,)7f + Cii^ . (4) 

The value of 77 in the vicinity of is: 

= a(Te - T)/{2C) . (5) 

From the expression for the thermodynamical potential all the interesting 
thermodynamical quantities and their jumps at the critical temperature can be 
derived. 

Coming back to the toy model, consider the total potential (eq. (||)) and 
expand it as function of 9 around the value = 0. The result is 

U = - {Rn sin af + R{g cos a- Rn'^ cos 2a) 0'^ 

+R[Rn^ COS 2a ~ {g cos a) /4] 9^ /i + 0{9^) . (6) 

It is easy to see the analogy between expressions (||) and (||). In fact 9 is the 
order parameter of the toy model: when 9 = the system is axisymmetric 
and it breaks this symmetry for 6* > 0. The potential U is the analogue 
of the thermodynamical potential <f>. No extra terms with odd power of 9 
appear, there is a critical value of for which the quadratic terms vanish 
and, last but not least, the quartic term is positive in the vicinity of ficrit- 
It is straightforward to define a "temperature" of the system. Let us put 
T = then Tc = cos2a /{gR cos a), a = 2(gi?cos a)^/ cos^ 2a and 

C = {gR cos a)^ / cos 2a. The reader can push the analogy further by defining 
the equivalent "thermodynamical quantities" and compute their discontinuities 
at the transition point. 

An equivalent treatment (but slightly more complicated) can be employed 
by replacing S7 by the angular momentum L. In this case a must be different 
from zero, in order to have a finite angular momentum when 9 = 0. 

2.2.2. The case of rotating stars 

The mechanism of axial symmetry breaking for a rotating star is very similar 
to the mechanism discribcd in the toy model. When the angular momentum L 
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of the star is larger than a critical value Lc, the symmetry of the star changes 
in order to increase its moment of inertia / and consequently to decrease its 
rotational kinetic energy for a given angular momentum L. 

In what follows, it is more convenient to parametrize the rotation of the 
star by the ratio R of the rotational kinetic energy to the potential energy: 
R = T/\W\ instead of the angular momentum L. In fact, the critical value 
of Rc at which the symmetry is broken depends only on the equation of state 
(hereafter EOS) of the fluid which constitutes the star. Moreover, in the case 
of a polytropic EOS, (P = an'^ , n and P being the density and the pressure 



respectively) Rc depends on 7 only. As recalled in § 2.1, for an incompressible 
star Rc = 0.1375 (secular instability). 

The symmetry breaking process is, of course, more complicated than the one 
described for the toy model. The number of degrees of freedom of the star is 
infinite instead of being just one for the toy model, consequently the analogy 
cannot be pushed further. 

In order to understand what follows, the following theorem is needed: 

Theorem: Among all the possible angular velocity distributions in a 
steady state configuration, rigid rotation (17 = const.) minimizes the 
energy of the star (see the problem at the end of this lecture) . 

We may now understand what happens to an axisymmetric rigidly rotating 
star when the angular velocity is slightly larger than the critical one. Consider 
for the moment the case without dissipative mechanisms. Since at constant 
angular momentum — larger than the critical one — the configuration which 
minimizes the total energy of the rigidly rotating star is a triaxial one, the 
star tends to break the axial symmetry and starts to deform to evolve to the 
new equilibrium configuration. In absence of viscosity or any other dissipative 
mechanisms, the angular momentum of each shell of the star is conserved and 
the angular velocity is no more uniform. It turns out that this new configu- 
ration, slightly non-axisymmetric has a total energy larger than the old one 
because the rotation is no more rigid (cf. the above theorem). Therefore the 
star does not evolve in this way and stays in the axisymmetric configuration. 

We may say that the conservation of the angular momentum of each shell 
acts as a potential barrier that inhibits the transition to the non-axisymmetric 
equilibrium configuration. The situation is completely different when viscosity 
is present. When the star starts to deform, the viscosity spreads out the angular 
momentum, so that the rotation becomes rigid and the star can go to the non- 
axisymmetric equilibrium configuration. It is obvious that the time taken to 
make the transition depends on the time required to rigidify the motion: the 
smaller the viscosity, the longer the time to make the transition. For this reason 
this instability is called secular. When the angular momentum is larger than a 
second critical value Lc,dyn the potential barrier due to the non-rigid rotation 
can be overcome and the transition toward the new equilibrium configuration 
is possible. The corresponding critical value of R is i?c,dyn ~ 0.2738 for an 
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incompressible fluid. This instability is called dynamical for it occurs on a 
dynamical timescale. 

The instability described above is at the origin of the transition of the star 
from a Maclaurin spheroid to a triaxial Jacobi ellipsoid. Since viscosity-driven 
evolution preserves the angular momentum but not the kinetic energy, the final 
equilibrium configuration has the same mass, the same angular momentum and 
rigid rotation. If the total mechanical energy is conserved but the conservation 
of the angular momentum is violated, then the Maclaurin steady state config- 
uration makes a transition toward the triaxial Dedekind ellipsoid, which is not 
rigidly rotating. 

Actually, there exists many other bifurcation points, all of them being parame- 
trized by a given value of the T/W ratio (pear configuration found by Poincare, 
fission bifurcation, etc..) 

From what has been said, the analogy between the mechanism of spontaneous 
symmetry breaking of the toy model and that of the rotating star should have 
clearly appeared. Consequently, the transition between the Maclaurin and 
the Jacobi ellipsoids can be studied by using the theory of second order phase 
transition, as we have done for the toy model. Historically, Bertin and Radicati 
were the first to recognize in 1976 that the bifurcations in the steady state 
configurations of rotating homogeneous stars are genuine second order phase 
transitions. More recently, Christodoulou, Kazanas, Shlosman and Tohline 
have systematically used this analogy to study the different bifurcations 
of rotating homogeneous stars and to investigate the possibility of a transition 
from one steady state configuration to another one. Not only have they found 
new results but they have also given more insight to the understanding of this 
mechanism. We highly recommend their articles to the 

interested reader. 



2.3. Previous results for compressible Newtonian stars 

The results on Maclaurin spheroids recalled above have been extended to com- 
pressible fluids, modeled by a polytropic EOS, by a number of authors. First 
of all Jeans , has shown that a bifurcation point towards triaxial con- 
figurations can exist only if the adiabatic index 7 is larger than 7crit — 2.2. 
The interpretation is that the EOS must be stiff enough for the bifurcation 
point to occur at an angular velocity lower than the maximum angular velocity 
r^K for which a stationary solution exists, fix is reached when the centrifugal 
force exactly balances the gravitational force at the equator of the star, and 
for this reason is called the Keplerian velocity; if the star were forced to ro- 
tate at > r^K, it would lose some matter from the equator. By numerical 
calculations, James [0 has refined Jeans' result to 

7crit = 2.238 (James 1964). (7) 

As concerns the secular instability in the compressible case, Ipser & Managan 
[2^ have shown that the m = 2 Jacobi-like bifurcation point has the same 
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location along uniformly rotating sequences as the m — 2 Dedekind-like point, 
as in the incompressible case. For 7 < 7critj the CFS instability still exists for 
modes I = m > 3 but not the viscous instability: if the viscosity is 

important, its effect is always stabilizing by acting against the CFS instability. 
Lindblom ||2^ has shown that for a star of 1.5 A/q constructed by means of 
a 7 = 2 polytropic EOS, the CFS instability is suppressed at temperatures 
T < 5 X 10® K by the shear viscosity and at T > 10^° K by the bulk viscosity. 
These results have been confirmed by Yoshida & Eriguchi (see also Sect. 2 
of ref. [0 and references therein). For completeness, let us note that if the 
neutron star interior is superfluid, the CFS instability is suppressed by the 
"mutual friction" of its components |^ . 

Regarding the gravitational wave signal from rotating neutron stars that un- 
dergo the above triaxial instabilities, Ipser & Managan |2^ have examined the 
case of polytropic stars with 7 = 2.66 and 7 = 3 which have bifurcated along a 
triaxial Jacobi-like sequence, under the effect of the viscosity driven instability. 
Wagoner has computed the gravitational signal from an accreting neutron 
star — modeled by nearly spherical homogeneous objects — for the five lowest 
modes of the CFS instability. Recently, Lai & Shapiro |^ have determined the 
gravitational wave form from newborn neutron stars — modeled as self-similar 
ellipsoids — undergoing the bar mode Z = to = 2 of the CFS instability. 

2.4. Generation of gravitational vi^aves 

It is worth discussing the role played by the two instability mechanisms, CFS 
and viscosity-driven, in the emission of gravitational waves by a rotating neu- 
tron star. Let us consider a just born neutron star (from the core collapse of 
a massive star) with a kinetic energy T larger than the critical one. For zero 
viscosity, the star breaks its symmetry via the CFS mechanism. Its fate is a 
triaxial Dedekind-like ellipsoid. Because the principal axes of this ellipsoid are 
fixed in an inertial frame the quadrupole moment is constant and the gravita- 
tional radiation stops. On the contrary if the viscosity is larger than a critical 
value , the CFS instability is inhibited and the star evolves toward a Jacobi 
configuration; it loses energy and angular momentum via gravitational radia- 
tion until its kinetic energy achieves the critical value Tc- It turns out that 
this second mechanism is more efficient for gravitational wave generation than 
the first one. In the intermediate case, the gravitational radiation light curve 
depends on the ratio between the rising time of the two instabilities. The real 
life is more complicated: in fact the EOS is far from being a polytropic. Con- 
sequently, the conditions for a bar instability depend on the EOS and on the 
mass of the star. Moreover, neutron stars are strongly relativistic objects and 
post-Newtonian or even post-post-Newtonian approximations are not sufficient 
to describe their steady states of rotation. A fortiori, the results obtained in 
the Newtonian theory can hardly be used to predict the gravitational radiation 
from these objects. 

In what follows, we shall try to answer to the following questions: 
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1. Among all the EOS proposed for nuclear matter, are there some stiff 
enough to allow the bar instability ? 

2. What is the influence of the general relativistic effects on the bar insta- 
bility ? 

3. In the case of a positive answer to the first question, and under the hy- 
pothesis that the general relativistic effects do not kill the bar instability, 
what are the masses of the stars for which the bar instability can develop ? 

4. Are these masses compatible with the masses of neutron stars already 
observed ? 



2.5. Finding the equilibrium configurations of a rotating star in the 
Newtonian regime 

In order to show how the Newtonian results can be extended to general relativ- 
ity, we need to explain how the equilibrium configurations of rotating barotropic 
stars are computed. Therefore, this section may appear quite technical. The 
reader not interested in mathematical and technical problems, can skip this 
section and go directly to section 2.9.2| . He should only remember that the 



exact solution for the motion of a triaxial rotating star in the framework of 
general relativity does not exist yet, the results presented in the next sections 
being obtained under adequate approximations. 

Consider a steady state star rotating with constant angular velocity fl about 
the z axis. In a frame rotating with the same angular velocity and about the 
same rotation axis, the star will appear static, therefore the problem of finding a 
steady state configuration becomes a problem of finding a static configuration. 
Let n, P{n) be the density and the pressure of the fluid respectively. The 
barotropic EOS P = P{n) is supposed to be known. The equilibrium equations 
reads 

l/ndP{n)/dp + dU/dp-n'^p = Q (8) 
l/ndP{n)/dz + dUldz^Q , (9) 

where U{x, y, z) is the self-gravitational potential, and = x"^ + y'^. By intro- 
ducing the specific enthalpy of the fluid, 

H{n) = f -^dn (10) 
J n dn 

the system of differential equations (^)-(||) can be integrated immediately to 
yield 

H + U - ]fl''^P^ = const. (11) 

Hence once U is known and H is fixed at the centre of the star, H and therefore 
n can be easily computed at all the points of the star. The surface of the star 
is defined by iJ = 0. 
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The potential U must be computed by solving in a consistent way the 2-D 
Poisson equation 

AU^inGn (12) 

where G is the gravitational constant. The integration of the system of equa- 
tions (12), ( [TTI ) and (10) is performed by relaxation: take a trial distribution of 



matter n, compute the potential U by solving the Poisson equation, compute 
the corresponding H by means of Eq. (|Tl|), inverse the relationship H{n) to 
obtain the density n and solve the Poisson equation to find a new potential 
U until convergence is achieved. About 50 iterations are necessary to obtain 
numerical solutions that differ less than lO"'^^ between two different iterations. 
In a such a way a 2-D equilibrium configuration is obtained. 

Once the equilibrium configuration is found, its stability can be studied. For 
this purpose we proceed in the following way: instead of stopping the relaxation 
at the J*'' iteration say, the procedure is continued altering slightly the gravita- 
tional potential at the ( J+1)"^ iteration: some perturbation e r^P2(cos 0) cos(2(^) 
is added to the axisymmetric potential, where = + z'^, P2 is the Legendre 
polynomial of degree 2 and e a numerical constant (e — 10^^). At the itera- 
tion J + 2 the perturbation is switched off and the relaxation continues. If the 
equilibrium configuration is stable, the perturbed solution relaxes to the ax- 
isymmetric solution; if not, it relaxes to the triaxial equilibrium configuration. 
In such a way, it is possible to find the critical value of 7, r/|iy|, etc... 

2.6. Extension to general relativity 

The analysis of the triaxial instability of rotating stars has been recently ex- 
tended to general relativity js^ in the case of rigid rotation. This latter as- 
sumption, crucial for finding a first integral of motion for non-axisymmetric rel- 
ativistic configurations by following the procedure of Carter jsj] , corresponds 
physically to the viscosity driven instability, but not to the CFS instability. 



The translation of the technique presented in § 2.5 from the Newtonian frame- 
work to the general relativistic framework is not completely straightforward. 
Before considering the relativistic analysis let us consider again the Newtonian 
description of triaxial configurations. 

2.6.1. Analysis of Newtonian triaxial configurations 

Let us first come back to the problem of finding axisymmetric stationary so- 
lutions of rigidly rotating stars in Newtonian theory. In this theory two ap- 
proaches are possible: The first one consists in finding an equilibrium configu- 
ration in a non- Galilean rotating frame (cf. § |2.5| ). In this frame the velocity 
vanishes and the term VP' p in Eq. (H) is the apparent acceleration (centrifugal 
acceleration). 

In the second approach (Galilean approach), the chosen reference frame is 
Galilean, the velocity of the fiuid f(0) = Q,p is not zero, and the term fip^ is 
due to the term v'^Vkv'^ of the Euler equation (inertial acceleration). In both 
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approaches the problem consists in finding a steady-state solution. Note that 
the mass conservation equation 

9n/at + V.(W) =0 (13) 

is identically satisfied and that the equation ( p^ ) for the gravitational potential 
U is the same one in the two approaches. 

The situation is different when looking for triaxial solutions. In fact for a 
Jacobi-like configuration, a steady state (equilibrium) solution exists only in the 
rotating frame. The questions that naturally arise are the following ones: is it 
possible to find an equivalent "equilibrium" configuration in a Galilean frame ? 
How to solve the Euler equation ? In order to answer these questions consider 
the unchanged Poisson equation (|l2|) and the Euler and mass conservation 
equations (x^ = p,x^ = z, = (f>) 

dv,/dt + v'^VkV^ + dH{p, z, (l},t)/dx' + dU{p, z,4>,t)/dx' = (14) 
dn{p,z,(l3,t)/dt + v''dn{p,z,(j),t)/dx'' + nVkv'' ^ (15) 

altogether with the Poisson equation ( [l^ ) for the potential which is unchanged. 
Let us look for a solution corresponding to a triaxial star undergoing rigid 
rotation about the z axis: U(p) = 0, ^(2) = 0, V(^^j = flp. The equations for the 
velocity read: 

-Q^p + dH{p,z,(j),t)/dp + dU{p,z,(l},t)/dp^O (16) 
dH{p,z,(l},t)/dz + dU{p,z,4>,t)/dz = (17) 
dH{p, z, (t>, t)/d4> + dU{p, z, (t>, t)/d4> = , (18) 

whereas Eq. ([l^) reduces to 

dn{p,z,(j),t)/dt + ndn{p,z,(f>,t)/d(j)^0 (19) 

It is easy to see that Eqs. (p^-(p^ are satisfied if H{p, z, cf), t) + U{p, z, cp, t) — 
1/2^^ p^ — C{t), where C{t) is an arbitrary function of t. This relation is a 
first integral of motion, very similar to relation (pT|). 

The continuity equation ( [l9| ) can be satisfied if we take C{t) = const and 
consider n, H, and U as functions of the new variables (p^z^ip) where 

i/j := (j) ~ fit . (20) 

Eq. ( p^ then becomes —fldn/dip + fldn/dip = 0, i.e. is identically satisfied. 

From a geometrical point of view, we may say that the Newtonian spacetime 
of a rotating triaxial star has a one-parameter symmetry group, whose trajec- 
tories (orbits) are defined by V = const. The spacetime vectors that generate 
this symmetry group are called Killing vectors. In the frame associated with a 
Killing vector, the star appears in a steady state. The use of the notion of sym- 
metry group or Killing vector fields allows the absolute definition (independent 
of any reference frame) of a steady state configuration. This is the definition 
that is going to be used in general relativity. 
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2.6.2. Rigid motion in general relativity 

Before the symmetry breaking, the spacetime generated by the rotating star 
can be considered as stationary and axisymmetric, which means that there 
exist two KiUing vector fields, fc" and m", such that is timehke (at least far 
from the star) and m" is spacelike and its orbits are closed curves. Moreover, 
in the case of rigid rotation, the spacetime is circular, which means that the 
2-spaces orthogonal to both fc" and are integrable in global 2-surfaces 
[ST) . This latter property considerably simplifies the study of rotating stars 
because some global coordinates {t, r, 9, 0) may be chosen so that the metric 
tensor components exhibit only one non- vanishing off-diagonal term {gt4>)- t 
and are coordinates associated with respectively the Killing vectors k" and 
m": k" = d/dt and m" = d/d(f>. The remaining coordinates {r,d) span the 
2-surfaces orthogonal to both k" and m". In all numerical work on rotating 
stars to date (see e.g. ref. for a review) isotropic coordinates are chosen, 
for which the two-dimensional line element differs from the flat space one by a 
conformal factor A^. In these coordinates, the components of the metric tensor 
are given by 



dx°'dx>^ = -N^ dt^ + B^r^ sin^ e{d(j3 - N'f' dtf + [dr^ + dO^] , (21) 



where the four functions iV, A^"^, A and B depend on the coordinates {r,9) 
only, the coordinates [t, (p) being associated with the Killing vector fields. 

When the axisymmetry of the star is broken, the stationarity of spacetime is 
also broken. As discussed above, in the Newtonian theory, there is no inertial 
(Galilean) frame in which a rotating triaxial object appears stationary, i.e. 
does not depend upon the time. It can be stationary only in a corotating 
frame, which is not inertial, so that the stationarity is broken in this sense. In 
the general relativistic case, where the notion of a global inertial frame is in 
general meaningless, a rotating triaxial system is not stationary for it radiates 
away gravitational waves. Even if a corotating frame could be defined, the 
body could not be in a steady state in this frame, because it loses energy and 
angular momentum via gravitational radiation. 

However, at the very point of the symmetry breaking, no gravitational wave 
has been emitted yet. As an approximation we neglect any subsequent grav- 
itational radiation. Then, for a rigid rotation, there exists one Killing vector 
field Z", which is proportional to the fluid 4- velocity u", hence 

M" = Ar, (22) 



where A is a strictly positive scalar function. Equation ( |22D is the definition of 
a rigid motion according to Carter W\ (i) there exists a Killing vector field; 



i^) More generally the relativistic concept of rigidity is defined by the requirement 
that the expansion tensor [of. Eq. below] associated with the vector field u" 



should vanish; it is easy to show that this requirement is fulfilled if takes the form 
(I22 
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(ii) the fluid 4- velocity is parallel to this Killing vector. In the axisymmetric and 
stationary case, the rigid motion corresponds to the constant angular velocity 
f2 := m'^/m*, the Killing vector entering equation ( p^ ) being then 

r = fc" + f7m", (23) 

where fc" and m" are the two Killing vectors defined above. The constancy 
of ri ensures that I" is a Killing vector too. The proportionality constant A 
of equation ( p^ ) is nothing else than the component u* of the 4- velocity u", 
where t is the coordinate associated with the Killing vector k". Note that the 
Killing vector I" is generally timelike close to the star and spacelike far from it 
(beyond the "light-cylinder"). 

In the non-axisymmetric case, k" and to" can no longer be defined as Killing 
vectors. We make instead the assumption that there exist (i) a vector field k" 
which is timelike at least far from the star, (ii) a vector field m", which com- 
mutes with k", is spacelike everywhere and whose field lines are closed curves, 

(iii) a constant fl such that the vector defined by the combination (^3|) is a 
Killing vector and (iv) the fluid 4- velocity u" is parallel to Z". These hypothe- 
ses are the geometric translation of the approximation of rigid rotation and 
negligible gravitational radiation. The commutativity of k" and m" ensures 
that a coordinate system {t,r,d,(j>) can be found such that k" = d/dt and 
to" — d/d(j). 



2.7. First integral of fluid motion in general relativity 

The purpose of this section is to show that in the framework of general relativity, 
there still exists a first integral for the fluid motion, this result being exact in the 
axisymmetric and stationary case, and being valid within the approximation 
of the existence of the Killing vector I" in the 3-D case. 

As detailed in ref. , the first integral can be found in an elegant geometri- 
cal way, using the canonical form of the equation of fluid dynamics introduced 
by Lichnerowicz js^ and Carter |Q and involving Cartan's calculus on differ- 
ential forms. However, in the present lecture, we follow a different approach, 
based on the formulation of relativistic hydrodynamics in an accelerated frame 
(Appendix In doing so, a parallel can be developed with the two points of 



view of Newtonian physics recalled in § 2.6.1 (that of a rotating observer and 
an inertial (Galilean) observer). 

The rotating observer is the observer Op comoving with the fluid, whose 
4- velocity is u" . The relativistic generalization of the "fixed" inertial observer 
is chosen to be the Eulerian observer Oe of the 3+1 formalism [Q, whose 
4- velocity is denoted by n". 

2.7.1. The point of view of the comoving observer 



Let us consider the Euler equation ( A21 ) in the frame of the fluid observer Of. 
Using notations of Appendix one has in this case, z;" = w", V" = and 
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E = e, so that Eq. (A21) reduces to 

^— V„j5 + a„ = , (24) 

e + p 

with tta — u^^f^Ua- Let us recaU that Vc denotes the covariant derivative 
within the 3-dimensional plane orthogonal to Of wordlines. Using Eq. ( ^ ) 
and the KiUing identity V(^alp) — 0, we get Oq — — VqIuA. Eq. ( |2^ ) then 
becomes 

^— V„p- V„lnA-0 , (25) 
e+p 

from which a first integral of motion is immediately obtained: 

i7 - InA = const. (26) 

where H is defined by 

H:= f^ = ln('-±P) , (27) 
J e+p \mBnJ 

niB being the mean baryon mass. The second equality in Eq. ( p7| ) is a con- 
sequence of the First Law of Thermodynamics at zero temperature. This def- 
inition of H constitutes a relativistic generalization of the specific enthalpy 
introduced by Eq. (|l^) . 

2.7.2. The point of view of the Eulerian observer 

Let us now consider the Eulerian observer Ce- His 4- velocity is 

n^^l-{k^+N^) , (28) 



where N is the so-called lapse function and A'^" the shift vector |38 . In the 
present case, N" = N'^ m" . Using notations of Appendix one has — n", 
i^ap = (the world lines of Oe are orthogonal to the hypersurfaces t ~ const) 
0ap = —Kgg, where Kap is the extrinsic curvature tensor of the hypersurfaces 
t = const. ||, V" = N-^ {n - N'l') m° = A^-^ - n", F = Nu* = NX. The 
Euler equation ( A19 ) then becomes 

1 



E + p 



V p + y"n^V^p a" = . (29) 



After straightforward calculations, making use of the fact that P is a Killing 
vector, the various terms that appear in this equation can be expressed as 

n^^V^V^ = -1 (n'^V^ In A^) P - ^Ky^ (30) 



(■^) the fact that the rotation 2-form of the Eulerian observer vanishes, explains 
why Oe is sometimes called the locally non-rotatmg observer |39|. 
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+2N-^ K'^^F + N-"^ K^.jn'' 71°' (31) 
a° = v"ln7V (32) 
K^^V^'V'' = N-^K^a^r = -r-^n^V^r + (r-2 - l)n''V^lniV (33) 

Accordingly, Eq. ( p9| ) becomes 

V a{H + \tlN -\tlT) + N^^n^'V ^,{H + \uN ~\TLT)la = . (34) 

Taking the scalar product with ^" and using the fact that is a symmetry 
generator, leads to n^V ^{H + h\N — InF) = 0. Introducing this latter relation 
in (^), we get 

H + In AT - In r = const. (35) 

This first integral of motion is exactly the same as that obtained in the co- 
moving frame [Eq. (p6|)] (remember that A = T/N). In the Newtonian limit, 
InN and InF ^ p^^^V2, so that Eq. (H) reduces to Eq. (|ll|), as expected. 



2.8. Gravitational field equations 

Let us now say a few words about the equations for the gravitational field 
in the general relativistic case. For axisymmetric rotating stars, these equa- 
tions are well known (see e.g. reL ||3^). As concerns the 3-D case and within 
the approximation made above (existence of the Killing vector P), a natu- 
ral approach would be to use the Geroch formalism pO| which reduces the 
4-dimensional Einstein equations to 3-dimensional equations by forming the 
quotient of spacetime by the trajectories of the Killing vector field. However 
in the present case, the Killing vector P is timelike inside the light cylinder 
and spacelike outside it. Consequently the second-order operators that appear 
in Geroch's formalism change from elliptic to hyperbolic type across the light 
cylinder. From the numerical point of view, we would rather have operators 
of a constant kind. For this reason, we do not consider Geroch formalism but 
instead the classical 3-1-1 formalism |3^, i.e. the foliation of spacetime by space- 
like hypersurfaces S f. Wi thin Sj, we choose coordinates (r, 61, 0) according to 
the prescription of § ^.6.2 . Next we introduce the coordinate ip as in Eq. ^ 



■0 := (j) — VLt. All the metric coefficients are then functions of the three variables 
(r, 9, ip). In the equations of the 3-f 1 formalism (cf. § 3.2 of |38)) appear partial 
derivatives with respect to t, r, 9 and (p. We replace them by partial derivatives 
with respect to r, 9 and ip according to the rule 

dt ' 91/- d(p ' dij' ^ ' 



We then obtain only 3-dimensional equations. Another simplification arises 
from considering only the dominant terms in the non-axisymmetric part of the 



16 



S. Bonazzola &: E. Gourgoulhon 



equations. We report to ref. p3| for further details. Let us simply mention 
that within our approximation, the triaxial metric tensor writes 

gaf3dx"dxf^ = -N{r,e,il;fdf + B{r,0,i^fr^sm'^e[d(j)-N'''{r,e)dtf 
+A{r, e, i^f [dr^ + dO^] . (37) 

Note that the metric coefficients given by Eq. ( p7[ ) are the components of the 
metric tensor with respect to the coordinates (i, r, 9, (p) and expressed as func- 
tions of the coordinates (r, 6, ip). 

Within this approximation, the Einstein equations reduce to four elliptic 
equations, for the functions A'', N'^ , AN and BN . In the Newtonian limit, the 
equation for N reduces to the Poisson equation 



2.9. Numerical results 



The non-linear 3-D elliptic equations resulting from the formalism developed 
above are solved iteratively by means of a spectral method The 
initial conditions are stationary axisymmetric configurations, constructed by a 
2-D general relativistic numerical code ||3^ , . The initial lapse function N is 
perturbed by a non-axisymmetric I ^ 2^m ~ ±2 "bar" term, in the same way 



as the Newtonian gravitational potential in § 2.5. The stability is determined 
by examining the subsequent evolution. 

2.9.1. Tests of the numerical code 

As a test of the relativistic code, the same value of the critical adiabatic index 
as that found by James |^ (Eq. |^) has been obtained at the Newtonian limit. 
Also in the Newtonian regime, it has been verified that when the polytropic 
index 7 increases, the ratio at which the triaxial instability occurs, 

tends to the clas sical value for incompressible fluids: r/|M^|crit(7 = 00) = 



0.1375 (cf. § 2A), as shown in Fig. g. Another test of the code consists in 
comparing results in the compressible polytropic case with previous numerical 
calculations, in the non-relativistic regime. Ipser & Managan ]43| and Hachisu 
& Eriguchi have obtained numerical models of Newtonian triaxial rotating 
polytropes, analogous to the Jacobi ellipsoids. They did not determine 7critlP] 
but performed calculations with fixed 7 > 2.66. The location of the triaxial 
bifurcation point along a sequence of 7 = 3 polytropes resulting from the 
relativistic code |Q used at the Newtonian limit has been compared with that 
obtained by the above authors. The agreement is better than 0.5% with the 
critical angular velocity of Hachisu & Eriguchi and of the order of 2% with 
their critical value of T/|iy|; with Ipser & Managan the agreement is 
better than 0.5% on both quantities. 



(*) However, Ipser & Managan state that their results indicate that the critical 
adiabatic index lies somewhere in the range 2.22 < 7crit < 2.28 
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Fig. 2. — Ratio of the kinetic energy T to the gravitational potential energy W at 
the triaxial Jacobi-like bifurcation point along a sequence of rotating Newtonian poly- 
tropes, as a function of the adiabatic index 7. The dashed horizontal line corresponds 
to the theoretical value of T/jM^I for incompressible Maclaurin spheroids. 



\ , I , \ 

0.8 0.6 0.4 



0.2 



Fig. 3. Critical polytropic index 7crit as function of the lapse function A'c measured 
at the centre of the star. Black boxes indicate configurations unstable with respect 
to radial oscillations. 
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2.9.2. Results for polytropes 

The investigation of rolativistic polytropic stars is interesting in several re- 
spects. First, it represents a natural extension of former classical works re- 
stricted to the Newtonian case. Second, a polytropic EOS does not suffer from 
the thermodynamical inconsistency relative to tabulated EOS (cf. Sect. 4.2 of 



ref. |42 ). It provides therefore an approximate but consistent model for real 
stars which allows a first investigation of relativistic effects. Since Newtonian 
polytropes obey a scaling law, 7crit5 the critical index for which the inset of 
the secular triaxial instability coincides with the maximum rotating case, is a 
global constant. In the relativistic case however, relativistic effects are sup- 
posed to influence the symmetry breaking and the critical index will depend on 
an appropriate parameter measuring the relativistic character of the object. A 
well suited quantity is the value Nc of the lapse function at the centre of the 
star. We thus obtain a 2-dimensional parameter space for maximum rotating 
stars which is intersected by a curve representing the metastable configurations 
and thus separating the regions of stable and unstable stars respectively. Fig- 
ure ^ shows the dependence of 7crit on Nc ranging from the Newtonian to the 
extreme relativistic regime. In the moderately relativistic domain there ap- 
pears a very slight decrease of 7crit- The negative slope reveals that the inset of 
relativistic effects tends to destabilize the star (within the limitations imposed 
by the approximate character of the theoretical approach [Q), the maximum 
decrease of 7crit being about 0.6%. In the strong field region we observe a 
smooth growth of 7crit beyond the maximum mass configurations due to the 
now persistently increasing stabilizing relativistic effects. 

2.9.3. Results for realistic equations of state 

We have determined in what condition the symmetry breaking may occur for 
rapidly rotating neutron stars built upon twelve EOS resulting from nuclear 
physics calculations. These "realistic" EOS are the same as those used in 
ref. and we refer to this paper for a description of each EOS. The EOS 
are labeled by the following abbreviations: PandN refers to the pure neutron 
EOS of Pandharipande [|5|, BJI to model IH of Bethe & Johnson |46|, FP 
to the EOS of Friedman & Pandharipande @, HKP to the no = 0.17 fm"^ 
model of Haensel et al. DiazII to model II of Diaz Alonso Glendl, 

Glend2 and Glend3 to respectively the case 1, 2, and 3 of Glendenning EOS 



|50 , WFFl, WFF2 and WFF3 to respectively the AVm + UVII, UV14 + UVII 



and UV14 + TNI models of Wiringa et al. (sT), and WGW to the AgJ,„„ + HV 
model of Weber et al. . 

Our results are shown in table ||. For a given EOS, the axisymmetric rotating 
models form a two parameter family; each model can be labeled by its central 
energy density Cc and its (constant) angular velocity f2. For a given value of 
Cc, n varies from zero to the Keplerian velocity flK- Following the method 
described above, we have searched for a symmetry breaking of configurations 
rotating at the Keplerian velocity. For five EOS (Glendl, Glend3, DiazII, BJI 
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EOS Afmax -^max -fbioak ^^c, break ^^^brcak 

[Mq] [Mq] [ms] [ms] [Mq] 



HKP 


2, 


.827 


3, 


.432 


0. 


.737 


1. 


.215 


0.161 


1.80 


WFF2 


2, 


T87 


2, 


.586 


0. 


.505 


0. 


.755 


0.30 


1.951 


WFFl 


2, 


T23 


2, 


.528 


0. 


.476 


0. 


.728 


0.27 


1.736 


WGW 


1, 


.967 


2, 


.358 


0. 


.676 




mar 


■ginally stable 


GlendS 


1, 


.964 


2, 


.308 


0. 


.710 






stable 




FP 


1, 


.960 


2, 


.314 


0. 


.508 


0. 


.604 


0.465 


1.736 


DiazII 


1, 


.928 


2, 


.256 


0. 


.673 






stable 




BJI 


1, 


.850 


2, 


.146 


0. 


.589 






stable 




WFF3 


1, 


.836 


2, 


.172 


0. 


.550 


0. 


.714 


0.325 


1.909 


Glendl 


1, 


.803 


2, 


.125 


0. 


.726 






stable 




Glend2 


1, 


.777 


2, 


.087 


0. 


.758 




mar 


■ginally stable 


PandN 


1, 


.657 


1, 


.928 


0. 


.489 






stable 





Table I. — Neutron star properties according to various EOS: M^lJ is the maximum 
mass for static configurations, M^^^ is the maximum mass for rotating stationary 
configurations, Pk is the corresponding Keplerian period, Pbreak is the rotation period 
below which the symmetry breaking occurs, -ffc, break is the central log-enthalpy at the 
bifurcation point and Mbrcak is the corresponding gravitational mass. The EOS are 
ordered by decreasing values of MjJai. 



and PandN), no symmetry breaking was found, whatever the value of JIk- 
These EOS are listed as "stable" in table |. For two EOS (WGW and Glend2) 
the evolution was not conclusive. These EOS are listed as "marginally stable" 
in table |. A better numerical precision could lead to a definitive conclusion. 
For five EOS (HKP, FP, WFFl, WWF2 and WFF3), the bar mode reveals 
to be instable for some Keplerian velocities. Table | gives the period Pbroak, 
gravitational mass Mbreak and central log-enthalpy ifc, break (cf. Eq. (|2^)) of the 
configuration having the lowest angular velocity and for which the symmetry 
breaking occurs. 

3. CW EMISSION FROM PULSARS 

As stated in the introduction, rapidly rotating neutron stars (pulsars) might 
be an important source of continuous gravitational waves. Morever, the ex- 
pected gravitational frequency is related to the rotation frequency and lies in 
the frequency bandwidth of the forthcoming LIGO and VIRGO interferometric 
detectors. While § |2| considers strong asymmetries resulting from an instability, 
we here investigate permanent slight asymmetries and the resulting continuous 
wave (CW) emission. 

Various kinds of pulsar asymmetries have been suggested in the literature: 
first the crust of a neutron star is solid, so that its shape may not necessarily be 
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Fig. 4. — Geometry of the distorted neutron star. 



axisymmetric under the effect of rotation, as it would be for a fluid: deviations 
from axisymmetry are supported by anisotropic stresses in the soUd. The shape 
of the crust does not depend only on the early history of the neutron star, 
especially on the phase of crystalization of the crust Q, but also on star quakes. 
Due to its violent formation (supernova) or due to its environment (accretion 
disk) , the rotation axis may not coincide with a principal axis of the neutron 
star moment of inertia and the star may precess. Even if it keeps a perfectly 
axisymmetric shape, a freely precessing body radiates gravitational waves. This 
efl^ect is discussed in P.Haensel's lecture [||. We consider another source of 
asymmetry here, linked to the magnetic field. Indeed, neutron stars arc known 
to have important magnetic fields and the magnetic pressure (Lorentz forces 
exerted on the conducting matter) can distort the star if the magnetic axis is 
not aligned with the rotation axis, which is widely supposed to occur in order 
to explain the pulsar phenomenon. 

3.1. Generation of gravitational waves by a rotating star 

Let us consider a neutron star rotating at the angular velocity f2 about some 
axis (cf. Fig. |^). In the framework of general relativity, the mass quadrupole 
moment relevant for gravitational radiation is Thome's quadrupole moment 
Xij Ipser has shown that the leading term in the gravitational ra- 

diation field hij is given by a formula which is structurally identical to the 
quadrupole formula for Newtonian sources (cf. § 2.2 of L. Blanchet's lecture 
jl)), the Newtonian quadrupole being simply replaced by Thome's quadrupole. 
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The non-axisymmetric deformation of neutron stars being very tiny, the total 
Thome's quadrupole can be hnearly decomposed into the sum of two pieces: 



7"rot I 7"! 



dist 



where is the quadrupole moment due to rotation {X\ 



rot _ 







if the configuration is static) and 1^^^^ is the quadrupole moment due to the pro- 
cess that distorts the star, for example an internal magnetic field or anisotropic 
stresses from the nuclear interactions. Let us make the assumption that the 
distorting process has a privileged direction, i.e. that two of the three eigen- 
values of Tf^^^^ are equal. Let then a be the angle between the rotation axis 
and the principal axis of 2"^"^* which corresponds to the eigenvalue which is 
not degenerate (cf. Fig. ||). The two modes /i+ and of the gravitational 
radiation field in a transverse traceless gauge are then given by ||55[| 



ho sin a 



hx = h 



sm a 



— cos a sin i cos i cos fl{t — to) 

1 + cos^ i 

— sma cos2r2(i — <o) 

— cosasini sinf2(t — to) 



■ sin a cos i sin 1^{t — to) 



(38) 



(39) 



where i is the inclination angle of the "line of sight" with respect to the rotation 
axis (see Fig. 0) and 

levr^G It 

= -^J^r ' ^^^^ 

where r is the distance of the star, P — 271/0, \s the rotation period of the star, 
I its moment of inertia with respect of the rotation axis and e := — 3/2Z|l''*// 
the ellipticity resulting from the distortion process. 

From the formulae (p8|)-(p^), it is clear that there is no gravitational emission 
if the distortion axis is aligned with the rotation axis (a = or tt). If both axes 
are perpendicular (a = 7''/2), the gravitational emission is monochromatic at 
twice the rotation frequency. In the general case (0 < \a\ < 7r/2), it contains 
two frequencies: Q and 20. For small values of a the emission at Q is domi- 
nant. Replacing the physical constants by their numerical values results in the 
following expression for the gravitational wave amplitude [Eq. (pO|)] 



ho = 4.21 X 10" 



24 


ms' 


2 


"kpc" 




\ ^ 1 




" e " 




[ P \ 




. r . 




Uo38 kgm^J 




[lO-eJ 



(41) 



Note that / = 10"^® kgm^ is a representative value for the moment of inertia of 
a 1.4 Mq neutron star (see Fig. 12 of ref. |56 ). 

For the Crab pulsar, P = 33 ms and r = 2 kpc, so that Eq. ( |4l| ) becomes 



7^ Crab 



1.89 X 10"^^ 



1038 kgm^ 



10- 



(42) 
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For the Vela pulsar, P — 89 ms and r — 0.5 kpc, hence 



h^""^^ = 1.06 X 10" 



27 


r ^ 1 








Uo38 kgm^J 




ilO-eJ 



(43) 



For the miUisecond pulsar^ PSR 1957+20, P = 1.61 ms and r = 1.5 kpc, 
hence 



1038 kgm^ 



10- 



(44) 



At first glance, PSR 1957+20 seems to be a much more favorable candidate 
than the Crab or Vela. However, in the above formula, e is in units of 10~^ and 
the very low value of the period derivative P of PSR 1957+20 implies that its 
e is at most 2 x 10~^ |5^. Hence the maximum amplitude one can expect 
for this pulsar is ~ 1.7 x lO'^^ and not 1.08 x IQ-^" as Eq. (||) might 

suggest. 



3.2. The specific case of magnetic field induced deformation 

Let us assume that the distortion of the star is due to its own magnetic field. 
It is then expected that the ellipticity e is proportional to the square of the 
magnetic dipole moment A4 of the neutron star 

. = /)£^At^ (45, 

In this formula, R is the circumferential equatorial radius of the star and (3 
is a dimensionless coefficient which measures the efficiency of this magnetic 
structure in distorting the star. In the following, we shall call (3 the magnetic 
distortion factor. Provided the magnetic field amplitude does not take (un- 
realistic) huge values (> 10"'^^ T), the formula ( psf ) is certainly true, even if 
the magnetic field structure is quite complicated, depending on the assumed 
electromagnetic properties of the fluid: normal conductor, superconductor, fer- 
romagnetic... 

Now the observed spin down of radio pulsars is very certainly due to the low 
frequency magnetic dipole radiation. M is then linked to the observed pulsar 
period P and period derivative P by (cf e.g. Eq. (6.10.26) of ref. |Q) 

. ,2 47r 3c3 IPP 
Mo sm a 

where a is the angle between the magnetic dipole moment A4 and the rotation 
axis. For highly relativistic configurations, the vector Ad is defined in the 



(*) We do not consider the "historical" millisecond pulsar PSR 1937+21 for it is 
more than twice farther away. 
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weak-field near zone (cf. Sect. 2.5 of ref. so is a. Inserting Eqs. 

(46) into Eq. (E^) leads to the gravitational wave amplitude 



6/3- 



crP sin a 

which can be cast in a numerically convenient form: 



6.48 X 10 



'30 



sin^ a 



\ ^ 1 


2 


"kpc" 




"ms" 




r P 1 


.10 km 




. r . 




P 




.10-13. 



im and 



(47) 



(48) 



Among the 706 pulsars of the catalog by Taylor et al. ||6^ , , the highest 
value of ho at fixed a, (3 and R, as given by Eq. (^^, is achieved by the 
Crab pulsar (P = 33 ms, P = 4.21 x 10-13, r = 2 kpc), followed by Vela 



(P = 89 ms, P = 1.25 x IQ-", 
P = 1.54 X 10-l^ r = 4.4 kpc): 



0.5 kpc) and PSR 1509-58 (P = 151 ms, 



F Crab 
'0 



'Vela 



1,1509-58 



I 1957+20 



4.08 X 10" 



= 1.81 X 10" 



1.50 X 10" 



31 


r R 1 


2 




(49) 




1.10 km J 




. 2 

sm a 


31 


r R 1 


2 


P 


(50) 




[ 10 km J 




SVC? a 


31 


r R 1 


2 


P 


(51) 




[ 10 km J 




siv? a 


37 


r R 1 


2 


P 


(52) 




1 10 km J 







We have added the millisecond pulsar PSR 1957 

-20 



3.1 



20 (P = 1.61 ms, P = 1.68 x 
to the list. From the above values, it 



lO"'^'^, r = 1.5 kpc) considered in 
appears that PSR 1957-1-20 is not a good candidate. This is not suprising since 
it has a small magnetic field (yielding a low P). Even for the Crab and Vela 
pulsars, which have a large P, the /iq values as given by Eqs. (|4^), ( ^ ) are, 
at first glance, not very encouraging. Let us recall that with the IQ-^^ Hz^^^^ 
expected sensitivity of the VIRGO experiment at a frequency of 30 Hz p^ , 
|63] , the minimal amplitude detectable within three years of integration is 



10" 



(53) 



Comparing this number with Eqs. (p9|)-(|50|), one realizes that in order to lead 
to a detectable signal, the angle a must be small and/or the distortion factor (3 
must be large. In the former case, the emission occurs mainly at the frequency 
f2. From Eq. (^) the gravitational wave amplitude can even be arbitrary large 
if a — > 0. However, if a is too small, let say a < 10"^, the simple magnetic 
braking formula ( ^ ) certainly breaks down. So one cannot rely on a tiny a to 
yield a detectable amplitude. The alternative solution is to have a large /3. For 
a Newtonian incompressible fluid with a uniform magnetic field, P — 1/5 p^ . 
In the following section, we give the P coefficients computed for more realistic 
models (compressible fluid, realistic equation of state, general relativity taken 
into account) with various magnetic fleld distributions. 



24 



S. Bonazzola &: E. Gourgoulhon 



Magnetic field lines 
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Fig. 5. — Magnetic field lines generated by a current distribution localized in the 
crust of the star. The thick line denotes the star's surface and the dashed line the 
internal limit of the electric current distribution. The distortion factor corresponding 
to this situation is /3 = 8.84. 



3.3. Numerical results 

We have developed a numerical code to compute the deformation of magnetized 
neutron stars within general relativity |55| . This code is an extension of that 
presented in ref. (59|. We report to this latter reference for details about the 
relativistic formulation of Maxwell equations and the technique to solve them. 
Let us simply recall here that the solutions obtained are fully relativistic and 
self-consistent, all the effects of the electromagnetic field on the star's equi- 
librium (Lorentz force, spacetime curvature generated by the electromagnetic 
stress-energy) being taken into account. The magnetic field is axisymmetric 
and poloidal. The numerical technique is based on a spectral method , . 

The reference (non-magnetized) configuration is taken to be a 1.4 Mq static 
neutron star built with the equation of state UV14 + TNI of Wiringa, Fiks & 
Fabrocini ]5l|| . This latter is a modern, medium stiff equation of state. The 
circumferential radius is i? = 10.92 km, the baryon mass 1.56 Mq, the moment 
of inertia / = 1.23 x 10'^* kgm^ and the central value of 500 is 0.36, which shows 
that such an object is highly relativistic. Various magnetic field configurations 
have been considered; the most representative of them are presented hereafter. 

Let us first consider the case of a perfectly conducting interior (normal mat- 
ter, non-superconducting). The simplest magnetic configuration compatible 
with the MHD equilibrium of the star (cf. ref |5^) results in electric currents 
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Magnetic field lines 
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Fig. 6. — Magnetic field lines generated by a current distribution exterior to a type 
I superconducting core. The thick line denotes the star's surface and the dashed line 
the external limit of the superconducting region. The distortion factor corresponding 
to this situation is /3 = 157. 



in the whole star with a maximum value at half the stellar radius in the equa- 
torial plane. The computed distortion factor is /3 = 1.01, which is above the 
1/5 value of the uniform magnetic field/incompressible fluid Newtonian model 
]5^ but still very low. 

Another situation corresponds to electric currents localized in the neutron 
star crust only. Figure |^ presents one such configuration: the electric current is 
limited to the zone r > — 0.9 rcq. The resulting distortion factor is /3 — 8.84. 

In the case of a superconducting interior, of type I, which means that all 
magnetic field has been expulsed from the superconducting region, the distor- 
tion factor somewhat increases. In the configuration depicted in Fig. ^, the 
neutron star interior is superconducting up to r, — 0.9 Teq. For r > r*, the 
matter is assumed to be a perfect conductor carrying some electric current. 
The resulting distortion factor is /? = 157. For = 0.95 r-cq, f3 is even higher: 
/3 = 517. 

The above values of /3, of the order 10^ — 10'^, though much higher than in 
the simple normal case, are still too low to lead to an amplitude detectable by 
the first generation of interferometric detectors in the case of the Crab or Vela 
pulsar [cf. Eqs. (|4^), ( pO[ ) and (|53|)]. It is clear that the more disordered the 
magnetic field the higher (3, the extreme situation being reached by a stochas- 
tic magnetic field: the total magnetic dipole moment M. almost vanishes, in 
agreement with the observed small value of P, whereas the mean value of 
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throughout the star is huge. Note that, according to Thompson & Duncan 



64 , turbulent dynamo amphfication driven by convection in the newly-born 
neutron star may generate small scale magnetic fields as strong as 3 x 10^^ T 
with low values of Bdipoie outside the star and hence a large /3. 

In order to mimic such a stochastic magnetic field, we have considered the 
case of counter-rotating electric currents. The resulting distortion factor can 
be as high as f3 — 5.7 x 10^. 

If the neutron star interior forms a type II superconductor, the magnetic field 
inside the star is organized in an array of quantized magnetic flux tubes, each 
tube containing a magnetic field Be ^ 10^^ T As discussed by Ruderman 
p5| , the crustal stresses induced by the pinning of the magnetic flux tubes is of 
the order BcB/2^q, where B is the mean value of the magnetic field in the crust 
{B ~ 10^ T for typical pulsars). This means that the crust is submitted to 
stresses ~ 10'^ higher than in the uniformly distributed magnetic field (compare 
BcB/2iio with i3^/2/io). The magnetic distortion factor /3 should increase in 
the same proportion. We have not done any numerical computation to confirm 
this but plan to study type II superconducting interiors in a future work. 



4. DISCUSSION AND CONCLUSION 
4.1. Spontaneous symmetry breaking 

From the results presented in Table |, it appears that only neutron stars whose 
mass is larger than 1.74 Mq meet the conditions of spontaneous symmetry 
breaking via the viscosity-driven instability. The above minimum mass is much 
lower than the maximum mass of a fast rotating neutron star for a stiff EOS 
(3.2 Mq [^). Note that the critical period at which the instability happens 
(P = 1.2 ms) is not far from the lowest observed one (1.56 ms). The ques- 
tion that naturally arises is : do these heavy neutron stars exist in nature ? 
Only observations can give the answer; in fact, the numerical modelling of a 
supernova core and its collapse ||66|| cannot yet provide us with a reliable an- 
swer. The masses of 17 neutron stars (all in binary systems) are known ||67[ |. 
Among them, four masses (all in binary radio pulsars) are known with a pre- 
cision better than 10% and they turn out to be around 1.4 Mq (see Table III 
in A. Wolszczan's lecture [Q). Among the X-ray binary neutron stars, two of 
them seem to have a higher mass: 4U 1700-37 and Vela X-1 (1.8 ± 0.5 Mq and 
1.8 ± 0.3 Mq respectively). These objects show that neutron stars in binary 
systems may have a mass larger than 1.7 Mq. 

A natural question that may arise is: why do X-ray binary neutron stars, 
which are believed to be the progenitors of binary radio pulsars, would have 
a mass larger than the latter ones ? We have not yet any reliable answer to 
this question. A first (pessimistic) answer is that the measurements of X-ray 
neutron star masses are not as reliable (compare the error bars of the masses 
of the binary radio pulsars with the ones of the X-ray binaries in Fig. 3 of 
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ref. |67|). Actually it should be noticed that the error bars of the X-ray pulsars 
do not have the same statistical meaning as the error bars of the binary radio 
pulsars p9| ] : they give only the extremum limits of neutron star masses in the 
X-ray binary. Consequently 1.4 Mq is not incompatible with these masses. 
If this is really the case forget all we have said about the above instability 
mechanisms: only the CFS mechanism for m > 2 can work, provided that the 
viscosity is low enough, which does not seem to be the case (especially if the 
"mutual friction" in the superfluid interior is taken into account [28| ). 

A related question arises naturally: why are the observed masses of millisec- 
ond radio pulsars almost identical ? Following the standard model, a millisec- 
ond radio pulsar is a recycled neutron star, spun up by the accretion of mass 
and angular momentum from a companion. The observed mass and angular ve- 
locity are those of the end of the accretion process. Consequently the accreted 
mass depends on the history of the system and on the nature of the compan- 
ion. By supposing "per absurdo" that all neutron stars are born with the same 
mass, it is difficult to understand why the accreted mass is the same for all 
neutron stars. A possible answer is that this could result from some obser- 
vational selection effect. For example, suppose that accreted matter quenches 
the magnetic field, it is then easy to imagine that the final external magnetic 
field depends on the mass of the accreted plasma. If the accreted mass is large 
enough, the magnetic field can be lower than the critical value for which the 
pulsar mechanism works. On the contrary, if the accreted mass is quite small, 
the magnetic field is large and the life time of the radio pulsar phase is shorter 
and consequently more difficult to observe. 

Let us conclude by saying that a lot of questions are still open about these 
systems. If we knew everything about neutron stars, such an observation would 
be a waste of time. The only thing that we can recommend is to stay open 
minded. 

4.2. CW emission from pulsars 

In this lecture, we have also investigated the CW emission resulting from the 
magnetic field induced distortion of neutron stars. The computations presented 
in § ^.3| show that the distortion at fixed magnetic dipole moment depends very 
sensitively on the magnetic configuration. The case of a perfect conductor inte- 
rior with toroidal electric currents is the less favorable one, even if the currents 
are concentrated in the crust. Stochastic magnetic fields (that we modeled by 
considering counter-rotating currents) enhance the deformation by several or- 
ders of magnitude and may lead to a detectable amplitude for a pulsar like the 
Crab. As concerns superconducting interiors — the most realistic configuration 
for neutron stars — we have studied type I superconductors numerically, with 
a simple magnetic structure outside the superconducting region. The distor- 
tion factor is then ~ 10^ to 10'^ higher than in the normal (perfect conductor) 
case, but still insufficient to lead to a positive detection by the first genera- 
tion of kilometric interferometric detectors. We have not studied in detail the 
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type II superconductor but have put forward some argument which makes it a 
promising candidate for gravitational wave detection. 

Acknowledgments 

We warmly thank Joachim Prieben and Pawel Haensel for their careful reading 
of the preliminary version of these notes. 

Problem: 

Show that, at fixed angular momentum, the kinetic energy of an incompress- 
ible fluid in a cylinder is minimal for rigid rotation. 

Solution: 

Let fl{p,z) be the angular velocity (p = \/x^ + y^). The kinetic energy of 
the fluid is given by T = 1/2 nVl'^{p, z)p'^ dV where n is the density of the 
fluid (assumed to be constant). We have to find an extremum of this quantity 
under the constraint that the angular momentum L = JynSl{p, z)p^ dV has a 
fixed value. By means of the Lagrangian multiplier technique, this amounts to 
find an extremum of 

/ \n^'^{p,z)p'^dV + \ [ nn{p,z)p'^dV , (54) 
Jv 2 Jv 

where A is the Lagrangian multiplier. By performing the variation with respect 
to Q. we obtain 

[ {n{p,z)p^ + Xp^)6ndv ==0 , (55) 

Jv 

from which O = const. 



Appendix A 

RELATIVISTIC HYDRODYNAMICS IN AN ACCELERATED 
FRAME 

In this appendix, we examine how to re-write the equation of momentum- 
energy conservation 

V^T"" = (Al) 

as a system of evolution equations with respect to a given observer O, the 
evolved variables being the energy density and the fluid velocity, both measured 
by O. We consider perfect fluids only: 

T"^ = {e+p)u"u^ +pg"^ , (A2) 

so that the equation for the fluid velocity will constitute a relativistic general- 
ization of the Euler equation of classical hydrodynamics. 
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The observer O is completely arbitrary; he is simply described by its 4- 
velocity w". Strictly speaking we consider a. family of observers O (a congruence 
of worldhnes), so that v°' constitutes a smooth vector field on spacetime. A 
fundamental tensor field related to O is the projection operator onto the 3-space 
P orthogonal to w": 

qa0 := 5a/3 + ""a W/3 ■ (A3) 

The 3-space P is made of spacelike vectors and can be thought as the "physical" 
three-dimensional space "felt" by the observer O. Note that if O is rotating 
7^ 0, see below), the vector space P is not integrable in global 3-surfaces. 
The motion of the observer O through spacetime is characterized by the 
Ehlers decomposition of V/jVq (see, e.g.. Sect. 4.1 of ref. [[70| ) 

VpVa = -t- Oap - aaVf3 , (A4) 

where 

:= qj^qp'^^iu Vf,] (A5) 

is the rotation 2-form of O, 

Oai3-^ qj'qp'^iuv^,) (A6) 
is the expansion tensor of O and 

fla -f^V^Wa (A7) 

is the ^-acceleration of O. 

For a 4- vector lying in the "physical space" P of O [v^j^W^ = 0), we 
introduce the 3-covariant derivative with respect to O as 

Vo^Wp := qJ^qp^'V^M. ■ (A8) 

This definition results in the following relation between the 4-covariant and the 
3-covariant derivatives: 

V^Wp = Va^Wp - vo. v^'V^Wp + {6e.^J.W^' - Lj^^,W^)vp , (A9) 

from which the following relation between the 4-divergence and the 3-divergence 
is immediately derived: 

V^W^" = V^M" + a^W' • (AlO) 
The fluid motion as seen by the observer O is specified by the Lorentz factor 

F := -v^ut" , (All) 

and the 3-velocity 

V" ■■= ^q'^X ■ (A12) 
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V" belongs to P and is the fluid velocity as measured by the observer O with 
his clock and his ruler. The following relations are immediate consequences of 
the above definitions: 



r = (1 - v^.v)-^/^ . 

The fluid energy density as measured by O is given by the formula 



or, according to the form (|A2|) of T^j^, 

E ^T^{e+p) -p 



(A13) 
(A14) 



(A15) 



(A16) 



From this expression, it is clear that the kinetic energy of the fluid with respect 
to O is included in the energy density E via the Lorentz factor F. 

Having set these definitions, let us now examine the equations of motion 
deduced from the momentum-energy conservation, Eq. (Al), which, using the 
perfect fluid form jA^ ) of Tq,^, can be written as 

(A17) 



(e + p)u''V,,u" + [(e + p)u^'] + V"p = 



The evolution equation for the fluid energy density E as measured by O is 
obtained by projecting Eq. (A17) along v". Invoking Eqs. (A13) and (All:), 
one obtains after straightforward calculations 

v'^V^E + V^ [iE+p)V^] + iE+p){2a^V^ + 6/ + O^^V^V) = . (A18) 

The relativistic generalization of the Euler equation is obtained by projecting 
Eq. ( A17 ) onto P, by means of qap. A fter s traightforward calculations (at a 
certain stage, use must be made of Eq. ( |A18| )) one obtains 



E 



(A19) 



Note that the flrst two terms on the left-hand side, w'^V^F" — a^y^w", con- 
stitute the Fermi- Walker derivative ]70| of V"" with respect to . The Fermi- 
Walker derivative measures the rate of change within P of V" with respect to 
the proper time of O. If the observer O has set up a local coordinate system, 
with respect to which the length of the vectors in P are evaluated, a derivative 
operator more convenient than w^V^ is the Lie derivative along v", £y. The 
term which naturally appears on the left-hand side of Eq. (A19) is then the 
convected derivative of V"' iTlfl, ||72| : 



:= £„ V 



{cu'^^+d")V^~a^V^v" . (A20) 
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The Euler equation (A19) then becomes 



^ ^V"p + D^p V") - a" . (A21) 



E + p 



In the Newtonian limit, reduces simply to d/dt. Moreover, the V^Vp,V°' 
gives the classical term {V ■ V)V and 2uj°'^V'^ gives the Coriolis term 2lj x V, 
induced by the rotation of the observer O with respect to some inertial frame. 
The terms involving 6af3 are due to the non-rigidity of the frame set up by the 
observer O. On the right-hand side, the classical [1/ p)Vp term is recognized, 
supplemented by the special relativistic term {dp/dt)V. The last term, the 
acceleration —a", contains gravitational as well as centrifugal forces. 
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